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SUMMARY 


PROBLEM 


Using the system performance equations developed in NUC TP 303, Handbook of 
Underwater Imaging System Design, develop a computer code for predicting the performance 
of underwater television systems. 

RESULTS 

A code called SCOUTS (for “Simulation Code for Underwater Television Systems”) 
and suitable for batch or demand runs has been written in FORTRAN IV. SCOUTS allows 
the user to select environmental, geometric, and hardware parameters. It can be used for 
system design as well as analysis. SCOUTS includes a parametric mode that allows automa- 
tic variation of multiple parameters. 
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1.0 INTRODUCTION 

This report documents the Naval Ocean Systems Center’s (NOSC) Simulation Code 
for Underwater Television Systems (SCOUTS). SCOUTS is a real time, interactive (when 
used in the DEMAND mode) computer code which is useful in the design, analysis and per- 
formance prediction of a large class of underwater television systems. A complete printout 
of SCOUTS in the FORTRAN IV language is included in the Appendix. For those with im- 
mediate applications, the information in Sections 2 and 5 is sufficient to describe the opera- 
tion and use of SCOUTS. Sections 3 and 4 give the theoretical basis for the calculations and 
internal program structure, respectively. The remainder of this Section concerns itself with 
the background and rationale of SCOUTS’ development. 

1.1 NEED FOR PERFORMANCE MODELING OF UNDERWATER TELEVISION 

During the past decade and a half man has rapidly extended his presence beneath the 
surface of the ocean. Manned submersibles, tethered vehicles, moored drill rigs and even 
underwater habitats reflect his desire to understand and utilize the unique world beneath the 
ocean’s surface. Underwater photography has been used throughout this era to document 
subsurface missions. This time period has also seen the development of underwater television 
to the point that it is equal in importance to photography as an undersea visual sensor. 

Modern underwater television systems use lenses, light sources and cameras specifi- 
cally designed for this environment. These components are incorporated into system con- 
cepts which attempt to minimize backscatter while maintaining an acceptable signal-to-noise 
ratio. As a result of the progress made during the last 1 5 years, the designer of underwater 
television systems now can choose from a wide variety of components and techniques in 
tailoring the system to his particular requirements. 

One rather disturbing effect of the availability of advanced components and tech- 
niques was that predicting the performance of the resulting system became much more 
difficult. Simple single scattering models gave erroneous results at the longer ranges corre- 
sponding to advanced imaging systems. More exact models were available, but they were 
primarily theoretical and not generally suitable for engineering use. Thus, by the early 
1970s, a situation arose where hardware was available which appeared to promise improved 
system performance, but the methodology for comparing the available hardware and accu- 
rately predicting performance was lacking. 

1.2 THE HANDBOOK OF UNDERWATER IMAGING SYSTEM DESIGN 

Under the sponsorship of the Navy’s Deep Ocean Technology Program, the Naval 
Ocean Systems Center undertook the task of providing a definitive handbook for the designer 
of underwater television systems. This effort produced the Handbook of Underwater Imag- 
ing System Design (HUISD) (Reference 1 ). HUISD provided a detailed discussion of optical 
water parameters, propagation of light underwater, characteristics of system components and 
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advantages and drawbacks of system concepts then in use or proposed. More importantly, 
HUISD gave a detailed design procedure for obtaining measures of system performance 
based on the above factors. This design procedure, validated by comparative system tests 
performed by NOSC at Morris Dam in 1970, was presented both as a set of nomograms and 
a step-by-step series of equations. These equations, which were more accurate and flexible 
than the nomograms, required lengthy hand calculations and were tedious to implement. 
Consequently, these system equations were not applied as widely as had been hoped, though 
their predictions proved accurate in those instances where the hand calculations were 
performed. 

1.3 EVOLUTION OF SCOUTS 

About a year after the publication of HUISD, NOSC began work on the computer 
modeling portions of the Advanced Unmanned Search System (AUSS) program. Part of the 
AUSS model required the simulation (i.e., performance prediction) of underwater television 
systems. Since the HUISD equations were already developed and validated, these equations 
were used as the basis for the AUSS computer code. Only the HUISD performance analysis 
of conventional systems (i.e., systems using non-laser sources and raster-scanned receivers) 
was coded, since these systems were the ones most commonly in use. The first AUSS TV 
program was coded in BASIC, but in order to achieve faster operating times the program 
was rewritten in the FORTRAN language. 

Experience with the AUSS TV program quickly showed that the HUISD system per- 
formance equations, coupled with real time, interactive demand terminal computing, pro- 
vided a powerful tool for TV system design and analysis. Using the program the operator 
was immediately furnished with the performance of the selected configuration as well as 
indications as to what (e.g., source power, backscatter) were the limiting factors. By succes- 
sively adjusting the TV system’s configuration the operator could optimize it for his particu- 
lar arrangement. This procedure was successfully used in designing and validating systems 
used in search, small-object detection, obstacle avoidance, etc. 

One of the problems with the AUSS TV program was that it is a subprogram integral 
to the much larger AUSS model. In order to use just the TV portion, much information ex- 
traneous to the operation of the TV subsystem has to be keyed in. Because there was a 
significant number of AUSS users who were interested primarily in the TV subprogram, it 
was decided to develop SCOUTS as an independent computer code for the simulation of 
underwater television. 

SCOUTS differs in a few important ways from the AUSS TV subprogram. SCOUTS 
is coded in FORTRAN IV and should be compatible as written with most medium- and 
large-size computers. A parametric capability allows several variables to be automatically 
varied according to limits defined by the user. Finally, although SCOUTS is primarily de- 
mand terminal-oriented, a batch option is included for computer facilities lacking real time 
capabilities. 



2.0 OVERVIEW OF SCOUTS 


This section will be concerned with defining the problem which SCOUTS solves, the 
• inputs required and the outputs provided. Coupled with the detailed operating instructions 

i in Section 5, this section provides sufficient information for successful operation of the 

program. 

1 2.1 SYSTEM GEOMETRY 

The geometry assumed by SCOUTS for its calculations is shown in Figure 1 . The 
source and receiver are spatially separated by the source-receiver separation, d. Without loss 
of generality the source is assumed to be to the receiver’s right when looking out towards the 
target. Their optical axes are coplanar and each is “canted in” by an angle 6 so that a perpen- 
I dicular dropped from the axes’ intersection point will bisect the line joining the source and 

receiver. The receiver’s (i.e.. TV’s) horizontal scan direction is parallel to the line joining 
source and receiver. 


The target’s center is in the plane containing the source and receiver’s optical axes 
(i.e.. the plane of the paper). The long dimension of the object is in the plane which is or- 
thogonal to the plane of the paper and parallel to the line joining source and receiver. The 




point of the object most distant from the source-receiver plane (i.e. the plane containing the 
source and receiver and perpendicular to the plane containing their optical axes) is assumed 
to be situated a distance R from it. The height of the object, H, is the extent of the object 
perpendicular to the source-receiver plane. 


The distance rg in Figure 1 is the distance past which backscatter contributes to the 
image. The depth of field, A, is the maximum distance from the target plane at which there 
is a “common volume” where source and receiver beams intersect any portion of the target 
whose height exceeds A will not be imaged. 


2.2 USER INPUTS 


Table 1 contains the inputs necessary to define the imaging problem for SCOUTS. 

All these variables must be entered on the initial calculation of any run. A rewrite capability 
(see Section 2.5) is included so that this information does not have to be entered each time it 
is desired to change only some of the variables. An example of the initial input for a demand 
run is shown in Figure 2. 

Table 1 . User inputs. 

Environment 

1 . Water type: coastal or deep 

2. Peak-to-trough wave height 

3 . Peak-to-trough bottom roughness 

Target Dimensions 

4. Target length 

5. Target height 

Platform Characteristics 

6. Source-receiver separation 

Source Characteristics 

7. Source type; incandescent, thallium iodide or mercury vapor 

8. Source input power 

9. Source full angle 

Receiver Characteristics 

10. Receiver type: vidiconorSlT 

1 1 . Receiver full angle 

12. Optics transmission 

13. F-number 
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CONVENTIONAL TV SYSTEM 


ENVIRONMENT 

COASTAL OR DEEP, 1 OR 2:? 

>1 

WAVE HEIGHT-PEAK TO TROUGH, FT:? 

>2 . 

BOTTOM ROUGHNESS-PEAK TO TROUGH, FT:? 

> . 5 

TARGET DIMENSIONS 
LENGTH OF TARGET, FT:? 

•10 . 

HEIGHT OF TARGET, FT:? 

>.5 

PLATFORM CHARACTERISTICS 
SOURCE RECEIVER SE PARAT I ON , FT : ? 

- 1 . 

SOURCE CHARACTERISTICS 

INCANDESCENT, THALLIUM IODIDE, OR MERCURY VAPOR; 1, 2 OR 3:? 
•2 

SOURCE POWER, WATTS : ? 

>200 . 

FULL SOURCE ANGLE, DEG:? 

•40 . 

RECEIVER CHARACTERISTICS 
VI DICON OR SIT, 1 OR 2 : ? 

>1 

FULL RECEIVER ANGLE, DEG:? 

•40 . 

TRANSMISSION OF OPTICS:? 

> .9 

f-number: ? 

>1.5 


Figure 2. Sample input for demand run. 


Tabic 2 is for the most part self-explanatory; liowevcr, a few comments are in order. 
SCOUTS was originally written for an object lying on the bottom viewed from a TV camera 
which was responding to ocean heave. To ensure an adequate depth of field. A, SCOUTS 
chooses A so that 

A = H-t-0.5 *(h, +h2) (1) 

where h j and h 2 are the wave height and bottom roughness respectively. For situations 
where the target is not on the bottom, hj and h 2 should be set to zero. 

SCOUTS incorporates information from HUISD to provide the appropriate character- 
istics for the water, lamp and receiver types. These characteristics will be discussed in detail 
in Section 3. However, in order for the user to supply meaningful inputs, it is necessary to 
discuss some of these characteristics here. 
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The two water types, coastal and deep, have optical characteristics identical with 
those plotted in Figures 6.2(B) and 6.2(C) of HUISD. The coastal water has a minimum 
attenuation, a, of 0.253/m (IM = 3.95 m) at X = 540 nm. For the deep water, the corre- 
sponding numbers are a = 0.049/m (1/a = 20.4 m) at 475 nm. The scattering coefficient, s, 
is 0.238/m for the coastal water and 0.030/m for deep water. 

For sources with conical beam patterns, the required “full source angle,” is just 
the apex angle of the in-water beam pattern. For a nonconical source, the “equivalent coni- 
cal beam pattern” (HUISD, p. 6-1 1, Eq. 6.10) must be determined and its apex angle used as 
input. For “source power,” SCOUTS requires the electrical power input to the lamp, then 
computes the light power output. 

The full receiver angle, 0 j^, is the full in-water horizontal field of view of the TV 
camera. SCOUTS assumes the standard 4-by-3 aspect ratio for the receiver’s field of view. 

For transmission of optics SCOUTS requires the overall decimal efficiency of the 
source and receiver. This is the combined efficiency of all the optical elements following the 
source and up to and including the camera lens. For example, if the collection efficiency of 
the light source and the transmissivity of the camera lens were each 0.9, the required optical 
transmissivity would be 0.81. 

The f-number, f/, is defined for small angles by 

f/4 (» 

where f is the focal length and D the receiver aperture diameter. The f/ has a minimum value 
of 0.5. Most camera lenses will indicate the f/ corresponding to each aperture stop. 

2.3 TYPES OF CALCULATIONS AVAILABLE 

2.3.1 Demand and Batch Runs 

SCOUTS can be run either as a batch or demand run. It is important to remember 
that SCOUTS was assigned for real time, interactive computing, i.e. for demand runs. De- 
sign problems are most efficiently handled in the demand mode, since the user will be able 
to refine his design based on immediate feedback of the performance of his previous choice. 
Another advantage is that after a period of running design problems via demand runs, the 
user will develop his intuition for good underwater systems design. 

The batch run capability was included primarily for computer systems not having 
real time facilities. However, the batch mode can be used to advantage in some situations. 
Some computer facilities charge substantially less for batch runs than for real time operation. 
Although SCOUTS’ costs are low because less than 1 .5 sec combined CAU (central arithmetic 
unit) and ER (executive request) time per calculation is required on ihe NOSC Univac 1 1 10, 
the cost savings might be significant if many calculations are required. Batch runs necr Barily 
provide a hard<opy output; this output may be absent or of poor quality on demand termi- 
nals In very large runs using the parametric mode (see Section 2.3.2) where many variables 
are to be parametized in a noninteractive fashion, using the batch mode allows the program 
to run withiHit the operator being present. 
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2.3.2 Nonparametric and Parametric Modes 

Nonparametric and parametric modes are available for both demand and batch runs. 

In the nonparametric mode, each input variable (see Table 1 ) takes on only one value and a 
calculation is immediately performed. In the parametric mode, the user selects one or more 
variables which he parametizes by assigning each a beginning value, incremental step and 
final value. SCOUTS then performs as many calculations as necessary to evaluate the system 
for all the selected values of the input variables, and the results are output. 

The parametric mode is thus similar to a succession of nonparametric calculations. 

The user selects the mode appropriate for his objectives. In optimizing a system, a succession 
of nonparametric calculations will allow the user to select new input values based on the re- 
sults of past calculations. The performance of a selected system for various values of input 
data can be documented more quickly in the parametric mode. 

2.3.3 Calculational Modes 

There are three calculational modes available to the user: single range, maximum 
range or maximum swath width. The calculational mode is chosen in demand runs by the 
appropriate response to a SCOUTS question and in batch runs by inclusion of the appropriate 
mnemonic (see Table 4). In the single-range mode the system performance is evaluated at the 
range supplied by the user. To understand the other two modes we must first define what is 
meant by usable swath width. 

The usable (or actual) swath width at range R is measured along that line in the 
source-receiver-target plane which is parallel to the line Joining source and receiver and a dis- 
tance R from it (the line AB in Figure 1 ). A segment of this line will have the following two 
properties: (1) a target located on this segment will liave an image contrast greater than or 
equal to 0.07 and (2) the target’s image will be large enough to span eight or more resolved 
TV lines. The length of this segment is taken to be the actual swath width. The actual swath 
width thus represents the lateral dimension in which successful imagery is possible. 

As a function of range the swath width increases almost linearly with range, reaches 
a maximum and then drops off rapidly as the range is increased further. The maximum range 
calculation evaluates the system at that range where the swath width, after having passed its 
maximum, drops to zero. Similarly, the maximum swath width calculation evaluates the 
system at the range which yields the greatest swath width. 

In the nonparametric mode. SCOUTS can evaluate the system in any of the three 
calculational modes. In the parametric mode, the user can have the system evaluated at 
either the maximum range or maximum swath width by indicating his choice and setting the 
range equal to zero. If the range has some other value the system will be evaluated at the 
selected range; however, SCOUTS will, at the user’s option, calculate the value of the maxi- 
mum range, maximum swath width, or both. Since SCOUTS performs these maximizing 
calculations by performing a series of single-range calculations, running time and costs can be 
minimized by selecting maximizing calculations only when necessary. 


II 



2.4 OUTPUT 


2.4.1 Minimum Range 

If the target is located too close to the source-receiver plane TV performance can be 
degraded. At distances which are too close, the resultant depth of field becomes less than 
the minimum value required by Eq. 1 . Additionally, at distances sufficiently short the source 
will be pointing directly towards the receiver and blind it. SCOUTS therefore calculates a 
minimum range prior to evaluating the user’s system. In demand runs this information is 
presented to the user in time for the user to choose a range greater than the minimum. 
SCOUTS then checks that the input range is indeed greater than the minimum. If not, the 
calculation will not be performed. 

One of the disadvantages of batch runs is that the user does not receive an indication 
of the allowable minimum range until the run is over. Thus, there is an uncertainty as to the 
smallest range that will execute sucessfully. Since SCOUTS sets the minimum range no 
smaller than the source-receiver separation, batch input ranges should be larger than this 
value. This will not guarantee a successful run but is the best the user can do with no other 
a priori information. 

2.4.2 Summary Output 

The results of a SCOUTS calculation (or calculations in the parametric mode) repre- 
sent the performance evaluation of the specified TV system. This performance evaluation is 
provided by SCOUTS in the form of two output summaries, the intermediate summary (Fig- 
ure 3) and the final summary (Figure 4). Intermediate summaries are provided at the end of 
each demand nonparametric run, while final summaries are output in all modes. Both sum- 
maries contain essentially the same information in different formats. The meaning of each 
of the outputs will now be presented in the order in which they appear on the final summary. 

The first thirteen columnar entires on Figure 3, from WATER TYPE to F-NUMBER, 
reproduce the user-selected input so that a record of the evaluated system is available along 
with the output. RANGE, FT and RANGE, AL give the range at which the calculation was 
performed in feet and attenuation lengths, respectively. All the performance measures from 
AVAILABLE LINES AT CTR to the second to the last entry refer to these measures evalu- 
ated at this range. MINIMUM, FT is the minimum allowable range (See Section 2.4.1 ) in 
feet. If the maximum range or maximum swath width options are chosen, the range at which 
each of these occurs is printed out in feet after MAXIMUM, FT and BEST, FT. respectively. 

The number of resolved TV lines spanning the object length is given by AVAILABLE 
LINES AT CTR for the object located at the center of the receiver’s field of view and by AT 
EDGE for the object located at the right-hand edge of the receiver’s field of view. Because 
of the interp-^lation subroutine used, the number of lines may be in excess of 525, which is 
taken to be the display limit. It has been shown that eight lines represent the threshold of 
object recognition (Reference 3). 

SCOUTS calculates the image contrast for an object having an inherent contrast of 
50 percent and a highlight reflectivity of 75 percent. CONTRAST AT CENTER and AT EDGE 
are respectively the image contrast for a target located at the center and right-hand edge of 
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CONVENTIONAL TV SYSTEM 


WATER TYPE: 



COASTAL 



SOURCE: THAL. lOD. 

WAVE HEIGHT, 

FT 

; 

2.00000 



SOURCE POWER: 

200.00000 

BOTTOM, FT: 



.50000 



SOURCE beam: 

40.00000 

LENGTH, FT : 



10.00000 



receiver: 

VIDICON 

HEIGHT, FT: 



.50000 



RECEIVER BEAM: 

40.00000 

TAU : 



.90000 






S.R. SEP.: 



1 . 00000 



F/ : 


1 . 50000 

RANGE, FT: 



24.40228 



RANGE, AL : 

1.87424 

AVAILABLE LI 

NES 

AT CTR. : 

393. 

90030 

AT 

EDGE : 

393.90030 

CONTRAST AT 

center: 

• 

16397 

AT 

EDGE : 

.07523 

AVAILABLE L 

H 

s 

SW WIDTH 

= 8. 

81963 

FT 

ACTUAL = 8 

.22299 FT 

LIMITED BY:R 

• 

ANG 






AVAILABLE R 

H 

S 

SW WIDTH 

= 8. 

95223 

FT 

ACTUAL = 8 

.22299 FT 

LIMITED BY:S 

• 

ANG 






AVAILABLE SW 

WIDTH = 17.7 

7186 

FT 


ACTUAL = 16 

.44598 FT 


Figures. Intermediate summary. 


the receiver’s field of view. Since backscatter is most severe towards the right edge of the 
field of view (remember that SCOUTS assumes the light is to the riglit of the receiver) the 
contrast at this point (as well as the number of lines available) will always be less than that in 
the center. 

The next eight entries on the final summary give information on the swath width to 
be expected. AV L.H.S. ( R.H.S.) SW WIDTH, FT is the available left-hand (right-hand) 
swath width. This is the segment, in feet, of the swath width line (See Section 2.3.3) to the 
left (right) of the center of the receiver's field of view which is illuminated by both source 
and receiver beams. ACTUAL, FT is the portion of the available left-hand (right-hand) swath 
width, in feet, where the target may be successfully imaged. Again, the three criteria used by 
SCOUTS for successful imaging are: ( 1 ) adequate depth of field as defined by Eq. ( 1 ), (2) at 
least eight resolved TV lines across the target’s longest dimension and (3) an image contrast 
of at least 7 percent. 

LIMITED BY gives a mnemonic wliich indicates what is limiting the extent of the 
R.H.S. (L.H.S.) actual swath width. The mnemonics S. ANG. and R. ANG refer to source 
angle and receiver angle respectively and indicate that the corresponding beam’s geometry is 
limiting the available swath width. BKSKTR indicates that the available swath exceeding 7- 
percent contrast is the limiting factor. Similarly, PWR signifies that the available swath width 
was limited by the requirement for eight TV lines/object length. DISPLY is printed when the 
number of lines/raster height required to place eight TV lines across the image’s longest 
dimension exceeds 525. 

AV SWATH WIDTH, FT is the arithmetic sum of the available left- and right-hand 
swath width and indicates the total available swath width for the given geometry and range. 
ACTUAL, FT is the corresponding sum of the actual swath widths, which is the linear 
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WATER TYPE: 

WAVE HEIGHT, FT: 

BOTTOM ROUGHNESS, FT: 
TARGET LENGTH, FT: 
TARGET HEIGHT, FT: 
S.R.SEP.,FT: 

SOURCE : 

SOURCE POWER, WATTS: 
SOURCE BEAM, DEG: 
RECEIVER: 

RECEIVER BEAM, DEG: 

TAU: 

F-NUMBER: 

RANGE, FT : 

RANGE, AL: 

MINIMUM, FT : 

MAXIMUM, FT : 

BEST, FT: 

AV LINES AT CTR; 

AT edge: 

CONTRAST AT CTR; 

AT EDGE: 

AV L.H.S. SW WIDTH, FT: 
ACTUAL, FT; 

- LIMITED by: 

AV R.H.S. SW WIDTH, FT: 
ACTUAL, FT: 

- LIMITED by: 

AV SWATH WIDTH, FT ; 
ACTUAL, FT ; 

BEST, FT : 


COASTAL 

2.0000 

. 5000 

10.0000 

.5000 

1.0000 

THAL. lOD. 

200.0000 

40.0000 
VIDICON 

40.0000 
. 9000 
1 . 5000 

24.4023 
1.8742 
2 . 5837 
. 0000 

24.4023 

393.9003 

393.9003 
.1640 
.0752 
8,8196 

8 .2230 
R. ANG 
8.9522 

8 .2230 
S . ANG 
17.7719 

16.4460 

16. 4460 


COASTAL 

2.0000 

. 5000 

10.0000 

. 5000 

1.0000 

THAL. lOD. 

200.0000 

40.0000 
VIDICON 

40 .0000 
.9000 

1 . 5000 
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Figure 4. Final summary. 


measure of the total horizontal distance over which satisfactory imagery is possible. BEST, 
FT is the maximum actual swath width. Zero will appear for this quantity unless the maxi- 
mum swath width option has been selected. 
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3.0 BASIS OF SCOUTS CALCULATIONS 

This section will outline the main features of the SCOUTS performance analysis and 
relate them to the corresponding steps in HUISID. This information is not strictly necessary 
for the general SCOUTS user, but it is useful to those who might desire to alter some of the 
program’s internal parameters to suit their specific problems. To that end, the various water 
and hardware characteristics used by SCOUTS will be reviewed in Sections 3.1 and 3.2. Sec- 
tion 3.3 will provide an outline of the calculational flow keyed to the labeled line numbers 
of the detailed listing in the Appendix. 

3.1 WATER CHARACTERISTICS 

All the water characteristics used by SCOUTS are documented and referenced in 
HUISID and appear as stored matrices entered via DATA statements in the program MAIN. 
The following discussion of the program’s features requires reference to the appendix, where 
the subroutines are arranged in alphabetical order. MAIN is the main program; the rest are 
subroutines. The function of each is explained in Section 4. 

The array A(l, J) (line MAIN 21) stores the attenuation coefficient for coastal (1=1) 
and deep (1=2) water for forward (J=l) and backscattered (J=2) light. It was obtained from 
HUISD (Table 6.1 2). The array B(1,J.K) (MAIN 24) contains the ratio of the effective to 
actual attenuation coefficient indexed according to water type (1), forward or reverse propa- 
gation (J) and full beam angle (K); these data also come from Table 6.1 2 in HUISD. The 
interpolation array of corresponding beam angles is stored in matrix T (MAIN 49). Scatter- 
ing coefficients, which were also obtained from Table 6.12, are stored in matrix S (MAIN 46). 
For both coastal and deep waters a rear hemisphere scattering (i.e., backscattering) percent- 
age of 2 percent is used (line TVFU 33). 

3.2 HARDWARE CHARACTERISTICS 

The spectral characteristics of the source, receiver and water enter into the HUISD 
systems performance analysis through the constructed functions r.(X) and H(X), which refer 
to forward and backscattered light, respectively. These functions are stored in the SCOUTS 
arrays G(1,J,K) and H(I.J.K) where 1, J, K refer to water, receiver and source types, respec- 
tively. These entires were computed by subtracting entries in Table 6.8 of HUISD in 
accordance with the effective bandwidths listed in Table 6.5. 

The imaging capability of the two types of cameras are stored in the array FM (I, J) 
(MAIN 30), where I refers to receiver type (i.e. 1=1= vidicon and 1 = 2 = SIT) and J indexes 
two constants per receiver type. For a given camera tube current, i, the number of resolved 
TV lines, N, is obtained according to 

gn (i/FM(l, 1 )) ( 3 ) 

FM(I,2) 
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The values of FM(I, J) were obtained by fitting Eq. (3) to available camera tube data. For 
the SIT the data were taken from HUISD, Figure 6.4(E), and the vidicon data from Reference 
2. In each case the 50-percent contrast curve was used. 

3.3 PERFORMANCE CALCULATIONS 

The first calculation SCOUTS performs for each system’s performance analysis is the 
evaluation of the minimum range. This occurs in lines HLIM 25-47. HLIMIT sets the mini- 
mum range equal to the greatest of (1 ) the source-receiver separation, (2) the minimum range 
which gives the required depth of field and (3) the minimum range which precludes the source 
and receiver looking at each other. 

All the other calculations are done in the subroutine TVFUNC and the routines 
called by it. After some initialization, TVFUNC sets the range equal to that desired for 
single-range calculation (line TVFU 53) or to 110 percent of the minimum range in the maxi- 
mizing modes (line TVFU 55). Next, a check is performed to make sure the display limit is 
satisfied (TVFU 59-64). If not, the remaining calculations are skipped and a flag is set so that 
a warning that the display limit has been exceeded is printed in the summary. 

In lines TVFU 72-120, SCOUTS evaluates the number of lines at the center and the 
edge of the field of view. If less than the minimum (i.e. eight) number of lines is obtained at 
the center of the field of view, SCOUTS assigns minimum left- and right-power-limited swath 
widths (TVFU 101-105). If there is no power limit even at the edge of the field of view, the 
maximum swath widths are associated with the power limit (TVFU 1 17-1 19). In the case 
where the eight-line limit is encountered at some angle other than the maximum or minimum, 
this angle is found in an iterative search (TVFU 87-99) and the appropriate swath widths are 
computed (TVFU 109-1 14). In every case, the photocathode current is obtained in subrou- 
tine DFUNC, which implements Eqs. 6.17 and 6.34 of HUISD. EFUNC evaluates the number 
of resolved TV lines, using the matrix FM(1, J). 

The effects of backscatter are accounted for in lines TVFU 121-1 27. Using the sub- 
routine XFUNC, which evaluates the photocathode current due to backscattered light (Eq. 
6.48 of HUISD), the contrast is calculated at the center and at the right-hand edge of the 
field of view. Using the interpolation routine, GFUNC, an appropriate right-hand swath 
width due to backscatter, W(7). is obtained. The geometric limitations on the swath width 
are next obtained (W(2), W(3), W(5), W(6)) through the use of UFUNC and VFUNC. FNK 
finds the limiting right and left swath widths and descriptors The limiting swath widths are 
summed to yield the total actual swath width. 

Lines TVFU 146-1 84 determine whether the calculation is for a single-range or a max- 
imizing calculation. If it is for a single-range calculation, TVFUNC returns contic! to HLIMIT. 
If not, TVFUNC determines whether the calculation is complete. If the last range increment 
is less than 5 percent of an attenuation length, control is returned to HLIMIT. If not, the 
range is incremented or decremented depending on the previous value of the swath width, 
and the performance is recalculated starting at line TVFU 40. 
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1 4.0 PROGRAM STRUCTURE 

4.1 SUBROUTINE ORGANIZATION 

SCOUTS consists of a main program and seventeen subroutines. The main program 
i defines constants used in performing the calculations and acts as a program driver. Eleven of 

( . the seventeen subroutines perform the computations. The remaining six subroutines - 

! SUBROUTINE DATAIN, SUBROUTINE CONVT, SUBROUTINE HLIMIT, SUBROUTINE 

' UPDATE, SUBROUTINE DATAOT, and SUBROUTINE SUMMARY were developed to 

[ provide flexibility of input and output. Figure 5 illustrates program hierarchy. A printout 

of each subroutine and of the main program is provided in the appendix. 

SUBROUTINE DATAIN reads user-supplied input data and sets program control 
Hags. If batch processing is desired, the user-supplied data are input via the card reader (unit 
number 5). Beginning in column one, the first card of this data deck must contain the five- 
‘ character alphanumeric word BATCH. This card sets the flag for batch processing. Subse- 

quent data cards are then read in an A6, 3X, 3F10.0 format. The alphanumeric information 
I is a mnemonic instruction that directs the program to a specific section of code. The data 

i values, where applicable, represent the initial, incremental, and final values of the data varia- 

I bles. Default values are set prior to reading any user-supplied data. 

1 

I 

i 



Figures. Program hierarchy. 
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If the program is being executed via a demand terminal, the terminal itself is the input 
device. In this case, interactive communication with the program is performed. The format 
of the input data, unless otherwise specified by the program, is GIO.O. It should be noted 
that the program can be executed via a demand terminal without interactive communication 
by typing the word BATCH as the initial response. Data are input as though batch processing 
were being performed. (For more detailed information on input data, see the section under 
program operation.) 

SUBROUTINE HLIMIT computes the minimum range and establishes the type of 
computations that will be made - single range, maximum range and/or maximum width. 

Incrementing of the input data variables is performed by SUBROUTINE UPDATE. 
This subroutine is called only when parametric variables have been defined. 

Two subroutines handle the printing of output data. Whenever execution is per- 
formed via a demand terminal, and no parametric variables have been defined, SUBROUTINE 
DATAOT prints the results after each evaluation is performed. This subroutine also writes a 
summary of the performance evaluation to the temporary disk file (unit number 9). At the 
conclusion of execution, SUBROUTINE SUMMARY reads the information stored on disk 
and prints tabular summaries of all performance evaluations computed. For batch processing, 
the printer (unit number 6) is the peripheral device utilized for output; for demand terminal 
processing, the terminal is the output device. 

SUBROUTINE CONVT was developed to allow the demand terminal user to input 
integer values left-justified. If a two-digit integer value has been entered, the second digit is 
read as alphnumertc data and “converted” to the appropriate numeric value. 

Data communication within the program is handled primarily through labeled com- 
mon blocks. The four blocks established - OPTS, TVCOM, lOLlST, and UPD — store user- 
defined control options, program computing constants, input/output variables, and input 
variable update information respectively. 

4.2 SUBROUTINE AND VARIABLE DEFINITIONS 

Tables 2 and 3 list and define the variables and subroutines used in SCOUTS. Table 2 
gives the type and definition for the variables that are listed in common statements and used 
by more than one subroutine. Although English units are used for input and output, MKS 
units are used within the subroutine calculations. Subroutine names, their functions and the 
programs called by them are shown in Table 3. 
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Table 2. Variable names and definitions. 


Descriptive Variable 

Type 

Definition 

NVIRON 

integer 

water type 

WAVEHT 

real 

wave height 

BTMRUF 

real 

buttom roughness 

TARGLN 

real 

target length 

TARGHT 

real 

target height 

HEIGHT 

real 

height (or range) 

SRSEP 

real 

source-receiver separation 

SPOWR 

real 

source power 

SDELF 

real 

full source angle 

IRTPE 

integer 

receiver type 

RDELF 

real 

full receiver angle 

OPTRAN 

real 

optics transmission 

FNUM 

real 

f-number 

ISTPH 

integer 

source type 

RANGMN * 

real 

minimum range 

RANGMX 

real 

maximum range 

RANGEB 

real 

best range 

SWWTHB 

real 

best swath width 

HTUS 

real 

usual height 

SWWTH 

1 

real 

actual swath width 

SWWTHL 

real 

actual L.H.S. swath width 

SWWTH R 

real 

actual R.H.S. swath width 

LIMFLG 

integer 

display limit flag 

NDEXWL 

integer 

L.H.S. width index 

NDEXWR 

integer 

R.H.S. width index 

LINCEN 

real 

available lines at center 

LINEDG 

real 

available lines at edge 

CTRCEN 

real 

contrast at center 

CTREDG 

real 

contrast at edge 


J 


i 



A 


: i 
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Tables. SCOUTS subroutines. 


Subprogram 

Function 

Reference 

Referenced by 


MAIN 

program driver; defines computing 
constants 

DATAIN, HLIMIT, 

DATAOT, UPDATE, 

SUMARY 

— 


DATAIN 

reads input data; sets control flags 

CONVT 

MAIN 


HLIMIT 

computes minimum range, defmes 
type of computation 

CFUNC, TVFUNC 

ATAN, TAN 

MAIN 


DATAOT 

prints results after each evaluation; 
writes summary data to disk 

- 

MAIN 


UPDATE 

increments parametric input variables 

- 

MAIN 


SUMARY 

prints tabular summaries 

- 

MAIN 


TVFUNC 

main computational subroutines 

DFUNC. EFUNC, VFUNC, 
XFUNC, CFUNC, UFUNC, 
FNK 

COS, ATAN, TAN 

HLIMIT 


CONVT 

converts a single character to 
numeric data 

- 

DATAIN 


XFUNC 

calculates backscatter current 

AFUNC 

TVFUNC 


AFUNC 

Simpson's rule integration routine 
for backscatter integral 

BFUNC, EXP. ALOG 

XFUNC 


BFUNC 

integrand for backscatter integral 

EXP 

AFUNC 

- 

CFUNC 

finds effective attenuation 
coefficient by interpolation 


HLIMIT 

i . 

DFUNC 

computes signal current 

TAN, ATAN, COS, EXP 

TVFUNC 

,, 

EFUNC 

evaluates number of lines 

ALOG 

TVFUNC 

i 

i 

FNK 

interpolates for limiting width 
and description 

- 

TVFUNC 

CFUNC 

interpolates for backscatter width 

TAN, ALOG 

TVFUNC 


UFUNC 

calculates l.h.s. receiver and r.h.s. 
source widths 

- 

TVFUNC 

r ■ ■ 

: i 

1 

VFUNC 

calculates r.h.s. receiver and l.h.s 
source widths 

- 

TVFUNC 

f i 
i I 
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5.0 PROGRAM OPERATION 


This section is intended to provide the user with information that is required for suc- 
cessful execution of the visual search system program. 

5.1 PRELIMINARY REQUIREMENTS 

The visual search system program utilizes three peripheral devices for transmission of 
data: the card reader, the printer, and a disk. The card reader and the printer are assigned 
the standard FORTRAN unit numbers 5 and 6 respectively. Unit 9 is assigned to a disk file 
used temporarily to store unformated data of summary output information. As required by 
the individual operating system involved, these unit numbers must be assigned to the appro- 
priate peripheral device prior to compilation of the program. 

After having assigned unit numbers to the appropriate peripheral device, the program 
is compiled and linked to create a set of absolute (or executable) binaries. Figure 6 illustrates 
the procedure used on the UNIVAC 1 1 10 to create the absolute element for execution. In 
the example the absolutes are stored under the name VSEARCH of the file TVSENS. The 
program driver and each subprogram are stored under separate file element names. 

@FOR,N TVSENS. MAIN 
@FOR,N TVSENS. DATAIN 
@FOR,N TVSENS. CONVT 
(3)FOR,N TVSENS. HLIMIT 
@FOR,N TVSENS. TVFUNC 
@FOR,N TVSENS. UPDATE 
@FOR,N TVSENS. DATAOT 
@FOR,N TVSENS. SUMARY 
@FOR,N TVSENS. AFUNC 
(3iFOR,N TVSENS. BFUNC 
@FOR,N TVSENS. CFUNC 
@FOR,N TVSENS. DFUNC 
@FOR,N TVSENS. EFUNC 
@FOR,N TVSENS. FNK 
@FOR,N TVSENS. GFUNC 
@FOR.N TVSENS. UFUNC 
@FOR,N TVSENS. VFUNC 
@FOR,N TVSENS. XFUNC 
@ PACK TVSENS. 

@PREP TVSENS. 

@MAP, N .TVSENS. VSEARCH 
IN TVSENS. 

END 


Figure 6. Creation of absolute element. 
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5.2 BATCH PROCESSING | 

i. 

(■ 

Batch processing of the absolute element VSEARCH is performed whenever the user 
does not wish interactive communication with the program. If batch processing is desired, 
the first card image of input data must be the Hollerith string BATCH. This five-character 
string informs the program that subsequent input will be in an A6, 3X, 3F10.0 format. The 
alphanumeric information is a mnemonic instruction that directs the program to a specific 
section of code. The data values, where applicable, represent the initial, the incremental, and 
the final values of the appropriate input variable. Table 4 lists the mnemonic instructions, 
the variable defined, and their default values. It should be noted that although the data as- 
signment instructions may occur in any order within the data deck, the STOP instruction 
must follow any given set of data assignment instructions. The last card image of the data 
deck must be the mnemonic instruction FINISH. Figure 7 illustrates a sample deck setup 
for execution of VSEARCH through batch processing. 

Table 4. Batch mnemonics. 



Mnemonic 

Defines 

Default Value 

NVIRON 

= 1., coastal; =2., deep 

1. 

ISTPE 

= 1., incandescent; =2., thallium; =3. mercury vapor 

1. 

IRTPE 

= 1., vidicon; =2., SIT 

1. 

WAVEHT 

wave height 

0. 

BTMRUF 

bottom roughness 

0. 

TARGLN 

target length 

0. 

TARGHT 

target height 

0. 

SRSEP 

source-receiver separation 

0. 

SPOWR 

source power 

0. 

SDELF 

full source angle 

0. 

RDELF 

full receiver angle 

0. 

OPTRAN 

optic transmission 

0. 

FNUM 

f-number 

.5 

HEIGHT 

range (or height) 

0. 

MAXSW 

sets flag to compute maximum width 

- 

MAXRAN 

sets flag to compute maximum range 

- 

SINGRN 

sets flag to compute single range 

- 

STOP 

stop reading input; end-of-record indicator for input 

- 

FINISH 

stop execution; end-of-file indicator for input 

j 

- 



BATCH 




NVI RON 

1 .00000 

.00000 

.00000 

WAVEHT 

2.00000 

.00000 

,00000 

BTMRUF 

.50000 

.00000 

.00000 

TARGLN 

10.00000 

,00000 

.00000 

TARGHT 

.50000 

. 00000 

.00000 

SRSEP 

1 . 00000 

.00000 

.00000 

1 STPE 

2 .00000 

.00000 

.00000 

SPOWR 

200.00000 

.00000 

.00000 

SDELF 

40 .00000 

.00000 

.00000 

IRTPE 

1 .00000 

.00000 

.00000 

RDELF 

40.00000 

.00000 

.00000 

OPTRAN 

.90000 

.00000 

.00000 

FNUM 

1 . 50000 

.00000 

.00000 

MAXSW 

.00000 

.00000 

.00000 

STOP 

.00000 

.00000 

.00000 

NVI RON 

2.00000 

.00000 

.00000 

RANGE 

65 .00000 

10.00000 

105.00000 

MAXSW 

.00000 

.00000 

,00000 

STOP 

.00000 

.00000 

.00000 

SRSEP 

1 . 00000 

1.00000 

4 .00000 

RANGE 

.00000 

.00000 

.00000 

MAXSW 

.00000 

.00000 

.00000 

STOP 

.00000 

.00000 

.00000 

SRSEP 

5 .00000 

.00000 

.00000 

SPOWR 

50 . 0 0000 

50.00000 

105 . 00000 

MAXSW 

.00000 

.00000 

.00000 

STOP 

.00000 

.00000 

.00000 

SDELF 

20.00000 

.00000 

.00000 

RDELF 

20.00000 

.00000 

.00000 

MAXSW 

.00000 

.00000 

.00000 

STOP 

.00000 

.00000 

.00000 

FINISH 

.00000 

.00000 

.00000 


Figure 7. Sample deck setup for batch execution of VSEARCH. 

5.3 DEMAND TERMINAL USE 

SCOUTS should be executed via a demand terminal whenever the user wishes inter- 
active communication. A demand run also provides intermediate information that can be 
useful in determining which input variables should be changed. 

Execution of the absolute element is performed by responding to the program’s 
questions. Questions concerning environment, source type and receiver type are answered in 
an 11 format. Other initial input is in a GlO.O format. After the initial input, SCOUTS asks 
for direction as to parametric or nonparametric, the calculational mode, whether another run 
is desired, etc. These questions are also answered in an 11 format. 

It should be noted that batch processing of SCOUTS can be performed from a de- 
mand terminal by responding ‘BATCH’ to the first question asked. In this case the succeed- 
ing input will have to follow exactly the format of Figure 7. 
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APPENDIX 


This appendix contains a complete FORTRAN IV listing of SCOUTS, arranged in 
alphabetical order by subroutine name. MAIN is the main program; all the rest are 
FORTRAN subroutines. All programs have compiled successfully on the NOSC UNIVAC 
1110 . 
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